home *** CD-ROM | disk | FTP | other *** search
/ Executor 2.0 / executorv2.0.iso / mac / demoware / Asymptote 1.1 Demo / Scripts / Demos / 9) Blackbody Figure / Plot_BlackBody_Figure < prev    next >
Encoding:
Text File  |  1994-11-01  |  1.7 KB  |  58 lines  |  [TEXT/DWat]
  1. % Plot a series of blackbody curves (Planck functions)
  2.  
  3. npts 200
  4. graphlimits 4 24 -30 15
  5. graphunits inches
  6. graphframe 1.5 7 2 9
  7. fullscreen yes
  8. graphticks -1 2 -1 2    % Log ticks, small every decade, big every 2 decades
  9. notation scientific    % Force axis labels to have only the exponent part
  10. digits 0        % in scientific notation
  11. fontsize 10        % Use a smaller font because there are so many labels
  12. drawframe
  13. drawxlabel Frequency, \sn [Hz]
  14. drawylabel Brightness B\d\sn(T) [erg s\\u-1\\d cm\\u-2\\d Hz\\u-1\\d steradian\\u-1\\d]
  15. scale scale 2 *        % Draw the title twice as big
  16. drawtitle Blackbody Curves
  17. scale scale 2 /
  18.  
  19. % Draw some double arrows to indicate basic wavelength bands
  20. scale scale 0.75 *
  21. drawdoublearrow xleft "ybottom 2 +" 11.5 "ybottom 2 +" Radio
  22. drawdoublearrow 11.5 "ybottom 2 +" 14.0 "ybottom 2 +" Infrared
  23. drawdoublearrow 14.0 "ybottom 2 +" 15.5 "ybottom 2 +" Optical
  24. drawdoublearrow 15.5 "ybottom 2 +" xright "ybottom 2 +" "X ray and ultraviolet"
  25. scale scale 0.75 /
  26.  
  27. % Restore default values
  28. fontsize
  29. notation
  30.  
  31. % Fill the x data with values that span the current graph limits
  32. xrange xleft xright
  33.  
  34. % Set a variable (T) and loop it to calculate each curve
  35. % The step expression "10 T * T -" increments T by a factor of
  36. % 10 with each step
  37. for T from 10 to 1E10 by 10 T * T - do Plot_Next_BlackBody T
  38.  
  39. % Add a second x axis to show the wavelengths
  40. notation scientific
  41. digits 0
  42.  
  43. fontsize 10
  44. % Offset the axis 3/4 of an inch below the first x axis
  45. graphframe gleft gright gbottom 0.75 - gtop
  46.  
  47. % Calculate the limits for the wavelength (lambda = c/nu) axis
  48. graphlimits c xleft 10^x / log c xright 10^x / log 0 0
  49.  
  50. % Draw the axis and its label
  51. drawaxis bottom
  52. drawxlabel Wavelength, \sl [cm]
  53.  
  54.  
  55.  
  56.  
  57.  
  58.